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Abstract 

Recent work that treats the Hawking radiation as a semi-classical tunnelling process 
from the four-dimensional Schwarzschild and Reissner-Nordstrom black holes is ex- 
tended to the case of higher dimensional Reissner-Nordstrom-de Sitter black holes. 
The result shows that the tunnelling rate is related to the change of Bekenstein- 
Hawking entropy and the exact radiant spectrum is no longer precisely thermal 
after considering the black hole background as dynamical and incorporating the 
self-gravitation effect of the emitted particles when the energy conservation and 
electric charge conservation are taken into account. 
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1 Introduction 



Hawking's remarkable discovery [1] in 1974 that black holes can radiate ther- 
mally reconciled a serious contradiction among General Relativity Quantum 
Mechanics and Thermodynamics at that time and put the first law of black 
hole thermodynamics on a solid fundament. At a big cost, however, this dis- 
covery also caused another controversial problem: what happen to information 
during the black hole evaporation? This notable problem is now often called 
as the information loss paradox [2], since the radiation with a pure thermal 
spectrum has no way to be recovered after black holes have disappeared com- 
pletely. Initially Steven Hawking and Kip Thorne thought that information is 
lost during the process of black hole evaporation, but John Preskill believed 
that information is not lost and can get out of the black hole. In 2004, Steven 
Hawking [2] changed his old opinion and argued that the information could 
come out if the outgoing radiation were not exactly thermal but had subtle 
corrections. 

Recently Parikh and Wilzcek [3] presented a greatly simplified model (based 
upon the previous developments by Per Kraus, et al. [4]) to implement the 
Hawking radiation as a semi-classical tunnelling process from the event horizon 
of the four-dimensional Schwarzschild and Reissner-Nordstrom black holes by 
treating the background geometries as dynamical and incorporating the self- 
gravitation correction of the radiation. The radiant spectra that they derived 
under the consideration of energy conservation give a leading-order correction 
to the emission rate arising from the loss of mass of black holes, which corre- 
sponds to the energy carried by the radiated quanta. Their result shows that 
the actual emission spectrum of black hole radiation deviates from strictly 
pure thermality, which might serve as a potential mechanism to resolve the 
information loss paradox. Following this program, a lot of recent work (see Ref. 
[5] and references therein) has been focused on extending this semi-classical 
tunnelling method to various cases of black holes such as those in de Sit- 
ter [6,7,8], anti-de Sitter [5,9,10] space-times, charged [8,11] and rotating [12] 
black holes as well as other cases [13]. But relatively less research [7,9,14] in 
higher dimensions is touched upon. 

In this Letter, we shall investigate the Hawking radiation via tunnelling from 
higher dimensional Reissner-Nordstrom black holes in de Sitter space-time. 
The subject is mainly motivated by the fact that recent brane- world scenarios 
[15] predict the emergence of a TeV-scale gravity in the higher-dimensional 
theories, thus open the possibility to test Hawking effect and to explore ex- 
tra dimensions by making tiny black holes in the high-energy colliders [16] 
(LHC) to be running in the next year. Within this context, recently there is 
much attention [17] being paid to examine the various properties, especially 
Hawking radiation of the higher-dimensional black holes. Our interest on this 
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topic is also due to the following two important aspects: (1) In the light of 
recent astronomical observation data for supernova of type la and the power 
spectrum of CMB fluctuations, it has been suggested that our universe is in 
a phase of accelerating expansion [18] and will asymptotically approach a de 
Sitter space. That is to say, the cosmological constant in our present universe 
might be a small but positive one. This realization has sparked a sense of ur- 
gency in resolving the quantum-gravitational mysteries of de Sitter space. (2) 
A conjectured de Sitter/conformal field theory (CFT) correspondence [19] de- 
fined in a manner analogous to the very successful AdS/CFT correspondence 
[20] has been proposed recently that there is a dual between quantum gravity 
on a de Sitter space and a Euclidean conformal field on a boundary of the de 
Sitter space. Much work has been done on finding a holographic description 
of de Sitter space and establishing an analogous duality. Although there has 
been considerable success along this line, the proposed duality is still marred 
by various ambiguities. Hence as an intermediate step, it may be appropriate 
to reinforce our understanding of de Sitter space at a semi-classical level by 
considering the de Sitter radiation as a semi-classical tunnelling process. 

On the grounds of all described above, it becomes obvious that the further 
study of black holes in higher dimensions is of great important. Of particu- 
lar interest is the case of Hawking radiation of the higher dimensional black 
holes. The aim of this Letter is to present a reasonable extension of the Kraus- 
Parikh-Wilczek's semi-classical tunnelling framework in the four-dimensional 
spherically symmetric space-times to the case of higher dimensional Reissner- 
Nordstrom black holes in de Sitter space-time. Moreover, we shall investi- 
gate the tunnelling radiation of the massless uncharged particles and massive 
charged particles across the black hole event horizon and the cosmological hori- 
zon, respectively. In addition, their possible phenomenological implications to 
the brane-world emission rate is briefly outlined. 

Our Letter is organized as follows. After introducing the metric of higher 
dimensional Reissner-Nordstrom-de Sitter black holes and their Painleve-type 
extensions in Section 2, we discuss Hawking radiation of massless uncharged 
particles and massive charged particles as a semi-classical tunnelling process 
from the black hole event horizon and the cosmological horizon in Sections 
3 and 4, respectively, where the tunnelling rate and emission spectrum are 
explicitly computed in each case. The Letter is ended up with a brief remark. 



2 Higher dimensional Reissner-Nordstrom-de Sitter black holes 
and their Painleve-type extensions 

The line element of (n + 2)-dimensional Reissner-Nordstrom-de Sitter black 
holes with a positive cosmological constant A = n(n+ 1)/(2/ 2 ) and the electro- 
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magnetic one- form potential are given by [21,22] 



ds 2 = -fdt\ + j~ x dr 2 + r 2 dQ 2 n 




'n ' 



(1) 

(2) 



where 



f = f(M,Q,r) = l 



+ 



00 n Q 2 



16n 



(3) 



2(n - l)V n r 2n ~ 2 P ' 



Here M and Q are the mass and electric charge of the black hole, respectively. 
I is the curvature radius of de Sitter space, V n denotes the volume of a unit 
n-sphere dVt 2 n . (Units G n+2 — c — h — 1 are adopted throughout this article.) 

When M = Q = 0, the solution (1) reduces to the pure de Sitter space with a 
cosmological horizon r = I which may be very large according to the existing 
knowledge of the cosmological constant. When the positive mass M increases 
with the charge Q decreased or fixed, a black hole horizon occurs and increases 
its size, while the cosmological horizon shrinks. For a special value of M and 
Q, the black hole horizon will touch the cosmological horizon. Besides, when 
the cosmological constant vanishes, the black hole has an outer/inner horizon 
located at 



For the general case in higher dimensions (n > 2), the horizons are determined 
by the equation f(M,Q,rh) = which is of (2n)-order so it in general has 
(2n)-roots. Typically for n > 2, there will be three positive (real) roots of 
f(M,Q,r h ), with the outermost root describing a cosmological horizon, and 
the remaining pair describing inner and outer black-hole horizons. Just as did 
in the four- dimensional case [8] , here we are only interested in two distinct real 
roots of them: the largest one is the cosmological horizon (CH) r c , the smaller 
one corresponds to the event horizon (EH) r + of the black hole. The explicit 
forms of these solutions are not needed for our discussions made here and are 
not illuminating though a detailed analysis [23] may be of some interest. 

Obviously, a key trick to apply the semi-classical tunnelling analysis is to intro- 
duce a coordinate system that is well-behaved at the horizons. Thus to remove 
the coordinate singularities at the horizons, one must perform a Painleve-type 
coordinate transformation in arbitrary dimensions 




M 2 - 



nQ 2 1 



(4) 




dt R = dt^f 



/ 



dr , 



(5) 
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under which the electric potential (2) retains its previous form modulo a 
gauge transformation, and the line element (1) is transformed to the following 
Painleve-like form 

ds 2 = -fdt 2 ± 2yjl - fdtdr + dr 2 + r 2 dVt 2 n , (6) 

where the plus (minus) sign denotes the space-time line element of the out- 
going (ingoing) particles across the EH and the CH, respectively. The new 
metric (6) is very advantageous for us to investigate the radiation of parti- 
cles tunnelling across the horizons and to do an explicit computation of the 
tunnelling probability. 

First, let us work with the metric in the new form (6) and obtain the ra- 
dial geodesies of the massless uncharged particles that follow the radial null 
geodesies, which is given by 



— = ±lT\jl - f ~ ±K h (r-r h ) , K h = K h (M,Q,r h ) = - — 



(7) 



r=r h 



where the upper (lower) sign denotes the radial null geodesies of the outgoing 
(incoming) radiation from the EH (CH), namely, = r+, r c , respectively. 

Next, it is known that different from the radial null geodesies of the massless 
uncharged particles, the trajectory followed by the charged massive particles 
is not light-like, it does not follow the radially light-like geodesies when it 
tunnels across the horizon. By treating the charged massive particle as a sort 
of de Broglie s-wave, its trajectory can be approximately determined as [8,11] 



where the plus (minus) sign represents the radial geodesies of the charged 
particles tunnelling across the EH (CH), respectively. 

In the following, we shall investigate Hawking radiation via tunnelling across 
the EH and CH, respectively, and calculate the tunnelling rate from each 
horizon. The overall picture of tunnelling radiation for the metric is very com- 
plicated, because radiation is both propagating inwards from the CH and out- 
wards from the EH. A formal, complete analysis must consider both of these 
effects, and there would undoubtedly be scattering taking place between the 
black hole horizon and cosmological horizons [7]. To simplify the discussion, 
we will consider the outgoing radiation from the EH, and ignore the incoming 
radiation from the CH when we deal with the black hole event horizon. While 
dealing with the CH case, we shall only consider the incoming radiation from 
the CH and ignore the outgoing radiation from the EH. This assumption is 
reasonable as long as the two horizons separate away very large from each 
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other, since the radius of the cosmological horizon is very large due to a very 
small cosmological constant, and the black hole event horizon considered here 
is relatively very small because the Hawking radiation can take an important 
effect only for tiny black hole typical of 1 ~ 10-Tev energy in the brane-world 
scenario. On the other hand, because Hawking radiation is a kind of quantum 
effect, it can be neglected and may not be observed for an astrophysical black 
hole with typical star mass about 1OM . 



3 Radiation of massless uncharged particles as tunnelling 



Before proceeding to discuss Hawking radiation of charged massive particles as 
a semi-classical tunnelling process, let us first consider the uncharged radiation 
outgoing from the EH. We adopt the picture of a pair of virtual particles 
spontaneously created just inside the horizon. The positive energy virtual 
particle can tunnel out and materialize as a real particle escaping classically 
to infinity, its negative energy partner is absorbed by the black hole, resulting 
in a decrease in the mass of the black hole, at the same time, the size of EH 
will reduce and the radius of CH will enlarge. Since the emitted particle can 
be treated as a shell of energy u, Eqs. (6) and (7) should be modified when 
the particle's self-gravitation is incorporated. Taking into account the energy 
conservation only, the mass parameter in these equations will be replaced with 
M — > M — uj when the particle of energy uo tunnels out of the EH. So the radial 
null geodesies of the uncharged particles tunnelling out from the EH is 

f = 1 - yfl-f(M-u,Q,r) ^k + (M- to, Q, r+) (r - r+) . (9) 

The imaginary part of the action of the uncharged particles can be expressed 
as 

ImS = Im / / dP' T dr = Im / / — d(M - J) , (10) 
Jn Jo J ri Jm r 

where we have used the Hamilton equation 

and denoted r\ and 77 the locations of the EH before and after the particle of 
energy u tunnels out, respectively. 

Substituting Eq. (9) into Eq. (10) and switching the order of integration, the 
imaginary part of the action can be easily evaluated by deforming the contour 
around the single pole r = r' + at the EH and reads 
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lmS = lm / ^=d(M - J) 

J m Jn l-^Jl- f(M-u',Q,r) 

i-M-u) ft 

-Jm MM-u/.g.r^-" 

27T -{r]-rfj=-lAS EH , (12) 



where Seh = Ki r +/4 = 47rr™/ (nu n ) is the Bekenstein-Hawking entropy of the 



EH. 



Since the geometrical optical limit and the "s-wave" approximation can be 
used here, using the semi-classical WKB method the tunnelling probability is 
found to be related to the imaginary part of the action via T ~ e~ 2Im5 , so the 
tunnelling rate at the EH is 

T ~ e" 2Im5 = e ^ (r /- r " } = e AS *« , (13) 



where ASeh = Seh(M — u>,Q) — Seh(M,Q) is the difference of black hole 
entropy after and before the particle emission. Obviously, the derived emission 
spectrum actually deviates from pure thermality, perfectly generalizing those 
obtained in Refs. [3,4,7]. 

Consider now the case of the massless uncharged radiation incoming from 
the CH, Eqs. (6) and (7) should also be modified when the particle's self- 
gravitation is incorporated. After taking into account the energy conservation, 
the mass parameter in these equations will be replaced with M — > M + u when 
the particle of energy uj tunnels into the CH, and the radial null geodesies of 
the uncharged particles tunnelling into the CH becomes 

r = -l + yJl-f(M + u,Q,r) « -k c (M + u, Q, r c ) (r - r c ) . (14) 
Applying the Hamilton equation 



the imaginary part of the action of the uncharged particles can be computed 
in a similar manner by deforming the contour around the single pole r = r' c 
at the CH 
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rrfc r — (M+Lj) rjrr- 

\mS = -\m / — d(M + uj') 



-m r 



-(M+u>) rr fc dr 

= Im / / . d(M + uj') 

Jn c l-j\-f{M + u>,Q,r) 



-M 

r- (M+u>) n 

J-M UM + u'Or') d{M + UJ,) 



k c (M + uj',Q,r' c 

~ 7 (>■'},. c) • (i6) 



17W. 



where r ic and rj c are the locations of the CH before and after the uncharged 
particle tunnels into the CH, respectively. In terms of the Bekenstein-Hawking 
entropy of the CH: Sch = V n r^/A = Airr™ / \nu n ) , the tunnelling rate at the 
CH is 

T ~ e ~ 2Im5 = e^r^-?c) = e ^ScH ? (17) 

where ASch = Sch(M + u,Q) — Sch(M, Q) is the change of entropy of the 
CH after and before the particle tunnels into the CH. Once again, the actual 
radiant spectrum is no longer precisely thermal. 



4 Radiation of massive charged particles via tunnelling 



Now, we proceed to discuss Hawking radiation of charged massive particles 
from the EH. Since the emitted particle can be treated as a shell of energy 
uj and charge q, Eqs. (6) and (8) should be modified when the particle's self- 
gravitation is incorporated. Taking into account the energy conservation and 
the charge conservation, the mass and charge parameters in these equations 
will be replaced with M — > M — uj and Q — > Q — q when the particle of 
energy uj and charge q tunnels out of the EH. So the outgoing radial geodesies 
of the charged massive particle tunnelling out from the EH and the non-zero 
component of electro-magnetic potential are, respectively, 



f(M-u,Q-q,r) 
2y/l-f(M-u,Q-q,r) 

Q-q 

(n - l)V n r n - 1 ' 



K + (M-u,Q-q,r + )(r-r + ), (18) 

(19) 



When we investigate the tunnelling process of a charged massive particle, the 
effect of the electro-magnetic field outside the black hole should be taken into 
consideration. So the matter-gravity system consists of the black hole and 
the outside electro-magnetic field whose Lagrangian function —(l/A)F flu F> lL ' 



8 



is described by the generalized coordinate = (A t , 0,0,0). As the general- 
ized coordinate A t is an ignorable one, to eliminate this degree of freedom 
completely, the imaginary part of the action should be written as 



-{Pr,P At ) 



ImS = Im / ! (L-P At A t )dt = Im P [ * " (f dP'-A t dP' Ai 
Ju v ' Jn U(o, o) v 



dr 



, (20) 



where r\ and 77 represent the locations of the EH before and after the particle 
of energy uj and charge q tunnels out, and (PA t , P r ) are two canonical momenta 
conjugate to the coordinates (A t ,r), respectively. 



Substituting the Hamilton's equations of motion 



dH 



r = 



A t = 



dP r 
dH 



(r;A t ,P At ) 



dP At 



(A t ;r,P r ) 



dH\ (T , Au p At) = d(M-u) , 
dH\ (At .^ Pr) = A t d(Q - q) , 



(21) 
(22) 



into Eq. (20), and switching the order of integration yield the imaginary part 
of the action 



ImS* = Im 



Im 



rr f r (M-u>, Q-q) 



l(M, Q) 



<n\f-J) ( \J „ , d(Q q') 



r 

Jn J(. 

(M- u ,Q-q) f r f 2J1- f(M -u',Q - q',r) 



{n - l)V n r r - 



dr 
f 



J(M, Q) Jn 



<(M,Q) Jn f(M -iu',Q-q',r) 
(n-l)V n r^ diQ - q) 



d{M - u') 



dr . 



(23) 



The above integral can be evaluated by deforming the contour around the 
single pole r = r' + at the EH. Doing the r integral first, we find 



ImS* : 



(M-u, Q-q) 



7Y 



(M,Q) k + (M -u',Q - q',r' + ) 



d{M - uo') 



(n - l)V n r. 



7^d(Q - q') 



(24) 



Using the area-entropy formulae of the EH: Seh = V n r^/A = inr™/ (nu n ), one 
can prove an identity which is essentially equivalent to the differential form of 
the first law of black hole thermodynamics 
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K+(M,Q,r + ) 



dM 



Q 



{n - l)V n r n + 



-dQ 



= —rl l dr + = -tiofi, 



-dS, 



(25) 
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By means of this identity, we can easily finish the integration and obtain 



J r , u n nw„. v ' 



1 ,(M-u, Q~q) 1 

= -- dS = --AS eh ■ 26 

2 J(M, Q) 2 

The tunnelling probability for the outgoing radiation at the EH is still repre- 
sented by Eq. (13) with ASeh = Seh(M — u,Q — q) — Seh(M, Q) being the 
difference of black hole entropy after and before the particle emission. Obvi- 
ously the derived emission spectrum actually deviates from pure thermality. 

Finally we turn to calculate the tunnelling rate at the CH. Different from 
the tunnelling behavior across the EH, the charged particle is found to tunnel 
into the CH. Ignoring the effect of the EH and taking into account the energy 
conservation and electric charge conservation, the mass and charge parameters 
in Eqs. (6) and (8) will be replaced with M — > M + u and Q — > Q + q, and 
generally speaking, the size of EH will expand and the radius of CH will 
shrink when the charged particle tunnels into the CH. At this moment, the 
radial geodesies of the charged massive particle tunnelling into the CH and 
the non-vanishing component of electro-magnetic potential are 

r = - /B±M±M « -« C (M + u,Q + q, r e ) (r - r c ) , (27) 

At = -{n Q l)V n r^- (28) 



It should be stressed that the first law of thermodynamics is satisfied both at 
the EH and CH, respectively, but one must note that the definitions of energy, 
temperature, and electric potential receive an opposite sign in both cases [21]. 
With this issue in mind, the calculation of the tunnelling rate at the CH is 
almost similar to those did in the EH case. In this way, one can get 



-(-M-u>,-Q-q) 

t-Q) 



Q + q' 



\mS = — Im / / 

Jr tc J(-M- 

.(-M-u-Q-q) rr fc 2 J I - f(M + U> , Q + q' , T 



Im 



r. 

J Tic 



(-M-Q) 

Q + q' , (n , , 

(n — l)V n r n 1 



."^ + - /) - (n-V>- '" + ^ 

d(M + u') 



dr 
r 



f(M + uj',Q + q',r) 
dr , 



(29) 



where r ic and rf c are the locations of the CH before and after the charged 
particle tunnels into the CH, respectively. Carrying out the integration, we 
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have 



ImS = 



,(-M-u, -Q- ? ) 7T T 

A-Af.-Q) « c (M + w',Q + g',r / c ) [ ( + " } 
Q + q' 



1 AM+lo, Q+q) I 

■7, / dS' = --ASch 

2 J(M, Q) 2 



2tt 



noj., 



(30) 



in which we have used an identity similar to Eq. (25) that holds at the CH. 
Thus the tunnelling rate is still given by Eq. (17) where AS C h = S C h(M + 
id, Q + q) — Sch(M, Q) is the change of entropy of the CH after and before 
the particle tunnels into the CH. Once again, the actual radiant spectrum is 
no longer precisely thermal and fits well to the universal relation presented in 
the standard semi-classical tunnelling formalism [3,4]. 



5 Concluding remarks 



In summary, we have investigated Hawking radiation of the massless un- 
charged particles and massive charged particles as a semi-classical tunnelling 
process from the black hole event horizon and the cosmological horizon of 
higher dimensional Reissner-Nordstrom black holes in de Sitter space-time, 
respectively. Our results show that the tunnelling rate is related to the change 
of Bekenstein-Hawking entropy and the exact radiant spectrum is no longer 
precisely thermal after considering the black hole background as dynamical 
and incorporating the self-gravitation effect of the emitted particles when the 
energy conservation and electric charge conservation are taken into account, 
thus perfectly extending the Kraus-Parikh-Wilczek's semi-classical tunnelling 
framework in the four-dimensional case to the higher dimensional spherically 
symmetric case. 

The leading-order, energy-independent term in the emission spectrum is the 
well-known pure thermal spectrum, but the energy-dependent terms when 
computed to any desired order in u, are indicative of a "greybody" factor in the 
emission spectrum. This means that a deviation from pure thermality probably 
have an important effect on the number of the brane localized particles emitted 
per unit time. This issue deserves further research in the future. 
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